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We report results from a computer simulation study on the rotational ratchet effect in systems 
of magnetic particles interacting via dipolar interactions. The ratchet effect consists of directed 
rotations of the particles in an oscillating magnetic field, which lacks a net rotating component. 
Our investigations are based on Brownian dynamics simulations of such many-particle systems. We 
investigate the influence of both, the random and deterministic contributions to the equations of 
motion on the ratchet effect. As a main result, we show that dipolar interactions can have an 
enhancing as well as a dampening effect on the ratchet behavior depending on the dipolar coupling 
strength of the system under consideration. The enhancement is shown to be caused by an increase 
in the effective field on a particle generated by neighboring magnetic particles, while the dampening 
is due to restricted rotational motion in the effective field. Moreover, we find a non-trivial influence 
of the short-range, repulsive interaction between the particles. 
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I. INTRODUCTION 

Thermal Brownian ratchets or, as they are sometimes 
called, Brownian motors, are devices that are able to ex- 
tract directional motion from Brownian random noise [l[ . 
In these out-of-equilibrium systems, it is possible to rec- 
tify the Brownian noise into directional motion. In ther- 
mal equilibrium such a phenomenon cannot exist as the 
second law of thermodynamics would be violated @, Q . 

Thermal ratchet effects have been known for a long 
time. Recently, however, they are again gaining attention 
due to their possible applications in biological transport 
d, Q and nanotechnology d, 0] ■ 

Most of the research on Brownian motors has been 
focused on directed translational motion. Exceptions are 
recent studies of the so-called rotational ratchet effect 
in ferrofluids, which has been investigated theoretically 
as well as experimentally (Tlj . 

Ferrofluids are suspensions of ferromagnetic colloidal 
particles (with diameters of about lOnm or larger) in a 
carrier fluid such as water or oil [IH, [l3[ . These systems 
can be driven out of equilibrium by, e.g., an oscillating 
magnetic field. The ratchet effect reported in [H, ||| con- 
sists of a noise-driven directed rotation of the particles, 
which are exposed to a field without a net rotating com- 
ponent. The rotations of the particles are associated to 
an effective torque, which is transferred to the solvent 
medium. This latter torque is of macroscopic size, mak- 
ing the fcrrofluid ratchet effect experimentally detectable 

m 

The theoretical investigations so far have been per- 
formed on the basis of the single-particle Langevin- and 
Fokker-Planck equations [H, [l(J. Interactions between 
the ferromagnetic colloids have mostly been neglected @ , 
the argument being that the concentration of magnetic 
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particles is extremely small in many ferrofluids (volume 
fraction » 1%). In more concentrated samples, however, 
one would expect the magnetic dipole-dipole interactions 
between the particles to become important. Indeed, a 
well known effect is the chain formation of the particles 
triggered by the anisotropy, particularly the head-to-tail 
preference, of the dipole interactions. There is, to our 
knowledge, only one theoretical study in which the im- 
pact of the dipolar interactions on the ratchet effect has 
been investigated [lfjj . This study approximates the in- 
teractions on a mean- field level, i.e., all the particles ex- 
perience a homogeneous effective field. 



In the present paper, we will investigate the impact 
of the true dipolar interactions on a particle level, i.e., 
by Brownian dynamics (BD) computer simulations. In 
this way we can not only capture the full anisotropy and 
range of the interaction (which is known to be crucial for 
self-organization processes in dipolar systems [T^ - tl6| ). 
but also the fact that the particles are mobile. 



This paper is organized as follows: In Sec. UH we 
present the model and the simulation methods used 
throughout this study. The next section deals with the 
rotational thermal ratchet effect in non-interacting sys- 
tems. Here, we will investigate the angular trajectories of 
the particles and the influence of the strength of the noise 
and the external field. In Sec. MI 51 we will then turn to 
systems, in which the particles interact via a short-range 
repulsive and a dipole-dipole potential. We will show 
that dipolar interactions can enhance as well as suppress 
the ratchet effect and we will analyze the mechanisms be- 
hind these effects. Further, we show that the short-range 
isotropic, repulsive potential has a significant influence 
on the ratchet behavior. The paper is then closed with a 
brief summary and conclusions. 
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II. MODEL AND SIMULATION METHODS 



where 



In this study, we consider a three-dimensional system 
of dipolar colloidal particles that are immersed in a sol- 
vent. Only the dipolar particles are handled explicitly. 
As a model we use a dipolar soft sphere (DSS) potential, 
which is comprised of a repulsive part U rcp and a point 
dipole-dipole interaction part C/ D : 

C/ DSS ( r<i , Mi, Mi ) = E^fo) + C/ D (ry , w , M ,) (1) 

In Eq. (JTJ, rjj is the vector between the positions of the 
particles i and j, its absolute value, and fii is the 
dipole moment of the ith particle. The dipolar and re- 
pulsive interaction potentials are given by 



^DSS 



J7 D (r y ,/ii,/ij) = 



and 



U lep (r) = £/ bb (r) - t/ bb (r c ) + (r, 



dt/ ss 

0— 5 — (7c), 
dr 



(2) 



(3) 



respectively. Here, [7 rep is the shifted soft sphere poten- 
tial, where 



U bb {r) = Ae 



(4) 



is the unshifted soft sphere (SS) potential for particles of 
radius a. 

We investigate the system by making use of Brown- 
ian dynamics (BD) simulations. These are based on the 
translational and rotational Langevin equations [TJI, [13 , 
which are integrated twice over a time interval that is 
larger than the inertial relaxation time and small com- 
pared to the time on which the configuration changes II? 
[20j . This procedure results in the equations [TJI [2l|, [22 



Ti(t + At) =r<(i) 



DyFiAt + JzOftetii (5) 



and 



e,(t + At) =ei(t) 
1 



DfTiAt x ei(t) + J2D$At$ x e<(t), (6) 



which form the basis of our BD simulations. Equations 
(jSJ and (j6|) correspond to solving the Langevin equations 
in the overdamped limit. In the equations above, = 
/tij//i is the orientation of particle i. The conservative 
forces and torques are given by 



rj-i rjiDSS _r_ rjiCXt 



(7) 

(8) 



= ~IH x^V w U DSS (r ij)Wlft ) (9) 



TCXt TjCXt 
j — /Xj x i$ 



(10) 



with an external field B cxt . In Eqs. (JSJl and ((6j), I? J and 
-Dq^ are the translational and rotational diffusion con- 
stants, which are given by 



k B T 
inr\o 



R _ k B T 

U <A — a • 



(11) 

(12) 



where is the viscosity of the solvent. The quantities 
and are Gaussian random variables that satisfy 



<£>=0, (£)=(), 



(13) 



= = <€TCj> = o. (w) 

Regarding the external field B ext , we use the same 
ansatz suggested previously in Refs. d, Specifically, 
the field has a constant component in :r-direction and an 
oscillating, yet asymmetric component in y-direction. A 
suitable ansatz is given by 



B ext (i) 



B y [cos(uJot) + sin(2uj t + S)]e y . (15) 



The important point is that this field involves only os- 
cillations, but no full rotations, irrespective of the phase 
shift S. Nevertheless, it turns out that the particles can 
perform directed full rotations. We note that the ansatz 
(|15|l is, however, by no means the only field with which 
a ratchet effect can be realized. In fact a multitude of 
different fields are suitable, if certain certain conditions 
are met: B x must be non- vanishing and there cannot be 
a At such that B y {t) = B y (-t + At) (see Ref. [9] for a 
detailed discussion of this issue). 

Note that the particles we consider here are immersed 
in a solvent. However, the solvent is only taken into 
account implicitly and on a single particle level, i.e., the 
random noise and the diffusion constants do not depend 
on the configuration of the particles. 

We consider N = 500 particles in our simulation box 
with periodic boundary conditions. The long-range dipo- 
lar interactions are taken into account by using the Ewald 
summation method [23| . 

For convenience, we make use of the following reduced 
units: Field strength B* = (cr 3 /e) 1 / 2 i? ct (a = x,y); 
dipole moment fi* = (ec 3 ) -1 / 2 ^; torque T* = T/e; 
time t* = tD'o / a 2 ; temperature T* = k B T/e; position 
r* = r/er. 

In addition, we will employ the parameter A = n* 2 /T* 
measuring the dipolar coupling strength relative to k B T. 
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III. RESULTS 

A. The thermal ratchet effect in a non-interacting 

system 

As a background for our investigation of the impact of 
dipolar interactions, we discuss in this section BD simu- 
lation results for the rotational thermal ratchet effect in 
systems of non-interacting particles. To this end, we ana- 
lyze the trajectories of the particles under the influence of 
the external field (fTSj) as well as the corresponding torque. 
Further, we investigate the dependence of the ratchet ef- 
fect on the strength of the external field versus that of 
the noise. We note that some of the points discussed 
in this section have already been investigated in Ref. @ 
via numerical integration of the respective Fokker-Planck 
equation. Our present BD results supplement these pre- 
vious theoretical results. 

The systems we consider in this section are character- 
ized by a temperature T* = 0.2. The particles are driven 
by a field of frequency Lua 2 /Dq =15, y-field component 
B* — 1, and various values of B*. The dipole moment 
is set to fi* = 1, such that the dipole-field coupling is 
(iB v /k B T = fi*B;/T* = 5. 

Since neither repulsive nor dipolar interactions are 
taken into account in this section, the density can be cho- 
sen arbitrarily. With these choices of the variables, our 
results are easily comparable to the ones from Ref. @ (for 
the precise relations between the dimensionless variables 
in our study and those in see the Appendix). 

To start with, we show in Fig. Q] the mean orientation 
of the particles, 



ICQ 




FIG. 1. Mean orientation of the particles during one rota- 
tional period of the field for non-interacting systems with 
B% =0.1 (black lines) and B* x — 0.3 (gray lines). The al- 
and y-components of S(t) are indicated by solid and dashed 
lines, respectively. The dotted line shows the field component 
in y-direction By Xt * . 



The second quantity cf> depicted in Fig. [5] corresponds 
to the averaged value of the angles fa of all the particles 



1 N 

Em*)- 

3=1 



(17) 



In contrast to fa this average angle <f> increases mono- 
tonically, since individual fluctuations are smeared out. 



s(*) = ]vE e *(*)' ( 16 ) 

i=l 

for the external field (H]) with B* = 0.1, 0.3 and 6 = 0. 
As can be seen, S x is essentially constant, while S y fol- 
lows (with a phase lag) the oscillating component of the 
external field indicated by the dotted line in Fig.[TJ Inter- 
estingly, while S x is increased for the field with B* = 0.3 
over the field with B* = 0.1, S y remains essentially un- 
changed. We will later see that this is of crucial impor- 
tance for understanding the impact of interactions. 

The behavior of the mean orientation seen in Fig.[T]ap- 
pears essentially deterministic. The actual ratchet effect 
is illustrated in Fig. [51 where we plot two angles fa and fa 
The former is the angle that an (arbitrary) particle i en- 
closes with the a;-axis. It first remains close to a multiple 
of 2tt (indicated by the horizontal lines) corresponding 
to the particle oscillating around the x-direction of the 
field. This behavior continues, until a noise-induced full 
rotation (i.e., a crossing of a horizontal line) occurs. One 
also sees that the forward rotation, i.e., an increase by 
2-7T, occurs more often than the corresponding backward 
rotation. This illustrates the directional character of the 
ratchet effect. 
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FIG. 2. The angular trajectories in terms of the polar angle 
4>i of an individual particle and the system-averaged angle (j> 
at B% — 0.3 and 5 — 0. 

Irrespective of these differences between fa and fa 
Fig- HI clearly demonstrates that there is a net rotational 
motion in one direction. This corresponds to the pres- 
ence of a net torque. We calculated the net torque as an 
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average of the time-dependent torque over one period of 
the field, i.e., 



where 



(T) = - / T(t)dt 



1 w 



with 



B, = -V Mj ^ C/ DSS (r y -, Ml , /x,) + B c 



(18) 



(19) 



(20) 



Numerical data for the net torque that the particles ex- 
perience for B* = 0, 0.1, and 0.3 over the phase difference 
S is presented in Fig.[3][cf. Eq. ([15])]. It is seen that (T*) 




FIG. 3. Averaged torque during one rotational period of the 
field as function of the phase difference S for B* = 0, 0.1, and 
0.3 (and B* = 1). 

is indeed non-zero and (at 5 = 0) positive, reflecting the 
net rotation of the particles to the right. Increasing the 
phase difference 5, the value of (T*) changes and even 
assumes negative values. This implies that the particles 
can perform forward as well as backward rotations de- 
pending on 5. 

We now consider in more detail the dependence of the 
net torque on the strength of the constant field contribu- 
tion B*. At B* = 0, the external field simply performs 
an oscillation into the y-direction. In that case, no net 
torque can be observed. With only one direction distin- 
guished by the field, directional rotational motion sim- 
ply cannot occur (2j. For non- vanishing x-components 
the magnitude of (f*) depends strongly on the value of 
B*. As illustrated by Fig. [3l increasing B* from 0.1 to 
0.3 results in considerably larger torques (T*). We ex- 
plain this increase as follows: At higher values of B*, the 
particles are much more aligned into the x-direction of 



the field, and thus in the plane of the field. The latter 
point is illustrated by Fig. [TJ S x is considerably larger 
for B* = 0.3 than for B* = 0.1 while S y remains es- 
sentially unchanged. In other words, higher values of 
B* ensure that the particles remain within the plane of 
the field without dampening the oscillations of the dipole 
moments in the y-direction. 

It is well established that ratchet effects, in general, de- 
pend strongly on the strength of the noise relative to the 
deterministic contributions to the equations of motion 
0. For the present system, this interplay is illustrated 
in Fig. @J where we plot the torque as a function of the 
dimensionless temperature T*. Inspecting Eq. ([5]), we 
see that the temperature T* influences the strength of 
the deterministic torque (due to the field) alone, if the 
diffusion constant is kept constant. In other words, T* 
is a measure for the aforementioned ratio of conservative 
torques to random noise. Small temperatures correspond 
to systems that are dominated by deterministic torques, 
while large temperatures correspond to noise-dominated 
systems. 

In Fig.|H we can see that the ratchet effect is strongest 
for finite temperatures in the range T* w 0.05 — 0.2. This 
means that the ratchet effect looses in strength for too 
small or too large noise contributions. If the temperature 
is too small, the field dominates the rotational motion 
of the particles, which are effectively unable to perform 
rotations against the field. At large temperatures, on the 
other hand, the noise dominates such that the influence 
of the field becomes insignificant. However, without the 
non- equilibrium influence of the external driving field, the 
ratchet effect cannot exist @, Q . 
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FIG. 4. The net torque as function of the temperature for 
B* = 0.1, 0.3. 

In Ref . Q , the behavior of a single dipole in the oscil- 
lating field (IToT) was investigated on the basis of a Fokker- 
Planck equation. Consistent with our results, the authors 
of Ref. [9] found a maximum in the strength of the net 
torque at finite values of the noise intensity. Moreover, 
for the particular choice B* = 0.3, our results for the 
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mean orientation and the net torque [see Figs. [T] and [3] 
are in quantitative agreement with those in Ref. [9] (see 
the Appendix for the relation between the dimensionless 
units). 



ture system, the maximum of (T*) is approximately 30% 
larger than the torque at A = 0. For T* = 0.2 the max- 
imum is even more pronounced: The torque is increased 
by nearly 40% compared to the non-interacting system. 
For higher values of A, (T*) decreases continuously. 



B. Influence of the particle interactions 

In a ferrofluid, the particles interact with each other 
via short-range repulsive as well as dipolar interactions. 
These interactions can be neglected in strongly diluted 
ferrofluids, they do, however, become important when 
the density of the dipolar particles becomes higher. 

As has been shown in previous studies, particle inter- 
actions can indeed have a profound influence on ratchet 
effects [H [l(| H3] : For instance, in a translational ratchet, 
they can reverse the direction of the effect or even give 
rise to it in the first place @, 25] . The latter is also true 
for the rotational ratchet effect. It was shown in Ref. 
that dipolar interactions treated on a simple mean-field 
level can induce effective particle rotations despite a van- 
ishing field component B x . 

In the following, we choose a density of per 3 = 0.2 cor- 
responding to a dipolar fluid of moderate packing fraction 
r\ = np<T 3 /6 ~ 0.105 [26|. This choice ensures that the 
dipolar interactions play a crucial role at higher coupling 
strengths. The frequency of the oscillating field is again 
set to u)G 2 1 Dq =15. Regarding the interaction parame- 
ters, we consider a range of values for the dimensionless 
dipole moment p* and different values of the dimension- 
less temperature T*. In this way we can explore both, 
impact of the dipolar interactions ©, and those of the 
repulsive interactions ([3]). 

Note that while we vary fi* (and, thus, A), we keep the 
products n*B* and n*B*, i.e., the dipole-field coupling, 
fixed. To indicate the used field strength, we therefore 
use the notation B+ =B*/fi*. 

In the lower temperature systems (T* — 0.2) consid- 
ered here, we use B+ = 1. With this choice the inter- 
action strength between dipoles and field remains as in 
the previous section. In the systems with T* = 1, wc 
use proportionally stronger external fields with B+ = 5. 
The relative strength of the external field compared to 
the Brownian noise is then equal to the one in the low 
temperature case. 

In Fig.[5]we present results for the z-component (T*) of 
the averaged torque (T*) for systems with different dipo- 
lar coupling strengths A = p* 2 /T* at B+ = 0.5, T* = 1 
and B+ = 0.1, T* = 0.2. At A = 0, the particles interact 
with each other via the soft-sphere interaction but not 
via the dipole-dipole interaction. Note that pure soft- 
sphere interactions do not affect particle rotations, and 
thus they should not influence the ratchet effect. There- 
fore the net torque found at A = for the T* = 0.2 
system equals the one shown in Fig. [3] for S = 0. 

Starting from the non-interacting system, the net 
torque plotted in Fig. [5] increases up to A ps 1.5 (T* = 1) 
or A « 2.5 (T* = 0.2), respectively. In the high tempera- 
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FIG. 5. The z-component (T*) of the averaged torque for 
different coupling strengths A at T* = 1, — 0.5 (solid 
line, right axis), and T* = 0.2, B+ = 0.1 (dashed line, left 



:E4 



0.05 
0.04 
0.03 
0.02 
0.01 
0.00 



tf \ 


— (T z ) 

- a> 


■ V 

•V 





0.08 
0.06 t) 

CM 



0.04 



0.02 



■0.00 







6 8 10 12 

A 



FIG. 6. The net torque (TJ), the mean traversed angle 
{4>), and the mean traversed angle (</>)2d calculated from the 
torque [see Eq. (|23|l ] are shown for T* = 1 and = 0.5. 

The behavior of (T*) is also reflected by the change in 
the (system-averaged) polar angle 4> 



(4>) = 4>{t + t) - 4>{t) 



(21) 



during one rotational period of the field. As shown in 
Fig- El (4>) an d (T*) behave almost identically, which can 
be understood by looking at the BD evolution equation 
([tj]). This equation corresponds to numerically integrat- 
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mg 



k R T 



(22) 



Here, = x e, is the angular velocity of particle i and 
d is a random Gaussian variable. Note that a differential 
equation for can be obtained from Eq. (f22j) by using 
the definition of u>i and taking the vector product with 
e^. Equation © then corresponds to integrating the re- 
sulting equation. Assuming the rotational motion of the 
particles to be restricted to the plane of the field (denoted 
by the subscript "2d" ) and neglecting the random noise, 
we find 



(ft 



2d 



r_Dl 



(T z ) 



2,1 



(23) 



from integrating both sides of Eq. (|22j) over one period 
of the field. This equation relates the traversed angle 
of the particles to the average torque. The dotted line 
in Fig. |6] demonstrates this relation. We calculated the 
traversed angle via Eq. (|2"3")1 from the torque component 
(T z ). It is seen that this equation slightly underestimates 
the observed value of (cf>) , which can be explained by the 
fact that Eq. (|23|) holds strictly only for particle rota- 
tions in the plane of the field. Additionally, we neglected 
the random noise, which is expected to introduce fur- 
ther deviations from the observed relation. Nonetheless, 
Eq. (|2"3"|) captures the general behavior of (T z ) as function 
of A quite well. 

We now discuss the origin of the maximum in Figs. [5] 
and |U The initial increase in strength of the ratchet 
effect for increasing values of A can be understood by 
considering the average effective field 



B c 



B c 



f 

M 



M 

i=i j=£i 



(24) 



felt by the particles, where M is the number of particles 
considered and 



B 



dip _ 3rjj (jCjj ■ fij ) 



(25) 



This effective field is depicted in Fig. [7] for particles with 
A = 1.44, T* = 1 and A = 2.45, T* = 0.2. These temper- 
atures and coupling strengths roughly correspond to the 
maxima in the net torque and the averaged traversed an- 
gle [see Fig. [S]- The plots in Fig. [7] show that, due to the 
dipolar interactions, the effective field components B X R * 
and By are increased as compared to the components of 
the external field. The fact that an enhancement of the 
(effective) field acting on the particles can support the 
ratchet effect, is already suggested by our results for a 
non-interacting system in Sec. IIII Al As shown in Figs. [3] 
and |4l an increase in B x alone or in both, B x and B y , 
can lead to larger values of the net torque. 

For coupling strengths higher than A » 1.5 (T* = 1) 
or A « 2.5 (T* = 0.2), respectively, the magnitude of the 
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FIG. 7. Mean values of the x- and y-components of the av- 
eraged effective field B cff * over one rotational period of the 
field at A = 2.45, T* = 0.2, and B+ = 0.1 (top) and A = 1.44, 
T* == 1, and = 0.5 (bottom). The solid lines represent 
the respective external fields, i.e., B£ =0.1, By = 1 (top); 
B+ = 0.5, B+ = 5 (bottom) 



ratchet effect (as measured by (T '*)) begins to decrease. 
We relate this behavior to the increase of the ratio of 
conservative torques induced by the dipolar interactions 
and the external field relative to the strength of the noise. 
However, in contrast to the (corresponding) decrease de- 
scribed in Sec. IIII Al [see Fig. H], a large contribution to 
the torque now stems from the dipole-dipole interaction 
and not from the particle-field interaction. This can be 
seen in Fig. [51 where we compare the functions S x (t) and 
S y (t) [cf. Eq. (HU)] for a system at A = 1.44 (T* = 1) 
and a more strongly coupled one at A = 9 (T* = 1). It is 
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FIG. 8. Mean orientation of the particles during one rota- 
tional period of the field for A = 1.44 (black lines) and A = 9 
(gray lines). The x- and y-components of S are indicated by 
the solid and dashed lines, respectively. 

Here, B+ = 0.5 and T* = 1. 
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seen that S x is significantly larger for the latter system 
(S x fa 0.58) than for the former one (S x ps 0.16). Recall- 
ing the discussion in Sec. IIII Al (non-interacting system), 
one would thus expect the ratchet effect at A — 9 to be 
even larger than at A = 1.44. However, the amplitude of 
S y is considerably smaller. Indeed, the maximum of \S y \ 
for the strongly coupled system is about 0.31, while it is 
about 0.72 for the one with A = 1.44. 

This means that the particles at A = 9 are much more 
aligned along the ^-direction (i.e., the constant part of 
the field) without closely following the oscillations in y- 
direction. In conclusion, the behavior seen at A = 9 is 
in stark contrast to what is shown in Fig. [T] for a non- 
interacting system. There, an increase in B* does not au- 
tomatically damp out the oscillations in the y-direction. 
Consequently, the ratchet effect is increased rather than 
damped. 

Finally, we note that, for relevant values of A, the 
relative increase in the net torque is larger for the low- 
temperature system than for the high-temperature one. 
Indeed, Fig. [5] shows that the value of T* influences the 
entire behavior of (T*) as a function of A. Therefore, 
not only the dipolar, but also the short-range repulsive 
interactions between the particles have an impact on the 
ratchet effect. 
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FIG. 9. (a) Pair correlation functions of systems at different 
temperatures (T* = 1, 0.2) but identical dipolar coupling 
strengths A. (b) The local field b(r) [see Eq. for the two 
different systems. 

The sensitivity against T* can be explained by the fact 
that the soft-sphere interactions affect the effective dis- 
tance between the dipolar particles. In Fig [51a), the 
radial distribution functions of two systems at identi- 
cal coupling strength A = 4 but different temperatures 
(T* = 1 and 0.2) are shown. Judging from the position 
of the main peak, two neighbors are typically closer to 
one another in the T* = 1 system than in the low tem- 
perature one. As a result the particles experience a con- 
siderably stronger effective field in the high-temperature 
system. This is illustrated by Fig. Hfb), where we plot 



the function 




where S r is the surface of a sphere of radius r. From a 
physical point of view, the function b(r) corresponds to 
the local field in a;-direction that is generated by neigh- 
boring dipolar particles with distance r* from the central 
one. From Fig. HJb), we see that the local field at short 
distances is significantly increased in the T* = 1 sys- 
tem as compared to the field in the T* = 0.2 system. 
Hence, as argued above, the rotational motion is much 
more restricted in the former system resulting in a less 
pronounced ratchet effect at a fixed coupling strength. 



C. Relation to self-assembly 

It is well established that strongly coupled dipolar par- 
ticles can self-assemble into a variety of structures includ- 
ing chains, networks, and sheets [HI Oil. Moreover, for 
dense systems of dipolar spheres, theory and computer 
simulations predict a phase with spontaneous long-range, 
parallel (i.e., ferromagnetic) order [13] ■ It is therefore an 
interesting question whether these phenomena have any 
relevance in the context of the rotational ratchet effect. 

The answer from our present BD simulations is essen- 
tially negative. Indeed, at the conditions where we found 
an increase in the ratchet effect (T* = 1, < A < 3.5; 
T* = 0.2, < A < 5.5) there is no global parallel order. 
Moreover, significant local ordering of the particles only 
occurs for dipolar coupling strengths A > 9 (T* = 1), 
which is outside of the range where we observe enhance- 
ment of the ratchet effect. The structures seen in such a 
highly coupled systems are illustrated by the simulation 
snapshot in Fig. 1101 Similar to ferrofluids subject to con- 
stant, homogeneous external fields, the oscillating field 
favors chain formation of the particles. Systems of lower 
coupling strength lack any such order. In particular, no 
local order can be observed at A = 1.44, T* = 1, i.e., 
where the ratchet effect is maximal. 

Another interesting aspect is the (possible) impact of 
the ferromagnetic phase transition occurring at higher 
densities. In Ref. [lfj, this question was investigated on a 
mean-field level, where the particles experience a (mean- 
field) torque of the form 

K N 

where K is the coupling strength. For sufficiently large 
K, one finds a spontaneous ferromagnetic ordering of 
the particles resulting in an effective non- vanishing field 
component in a;-direction fioj j . Due to this net field, the 
ratchet effect can occur even in the absence of an external 
x-component of the field. 
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FIG. 10. Snapshot of a system at A = 9, T* = 1, and = 
0.5. 



We have searched for a similar phenomenon in our 
many-particle system. However, at the parameters con- 
sidered, we were not able to find a net particle rotation. 
Not even at high densities, where the dipolar soft spheres 
undergo a ferromagnetic phase transition [27| , did we de- 
tect such a rotation. 

This could be due to several reasons: First, the true 
effective field within the ferromagnetic phase is inhomo- 
geneous and typically much weaker than any average 
"mean-field" (this is also the reason that the isotropic- 
ferromagnetic transition in a true dipolar system occurs 
at much larger coupling strengths than those predicted 
by mean-field theory [la [27l] ) . 

Second, in a dense dipolar system, the orientations of 
the particles are strongly coupled over large distances. 
In other words, the dipole orientations are severely re- 
stricted, which further suppresses the ratchet effect. 



IV. CONCLUSIONS 

In this study, we have investigated the rotational ther- 
mal ratchet effect for non-interacting particles as well as 
particles interacting via long-range dipolar interactions. 

With our particle based simulations, we looked at the 
angular trajectories of the dipolar particles, which con- 
clusively illustrate the net rotating behavior of the driven 
particles. For non-interacting particles, we found that 
a finite ratio of deterministic torques to random noise 
yields a maximally pronounced ratchet effect. 

The main focus of this study, however, was the inves- 
tigation of the influence of dipolar interactions on the 
rotational behavior of the particles. In particular, we 
showed that dipolar interactions can have an enhancing 
as well as a dampening effect depending on the dipolar 
coupling strength A. The enhancement found at small 



values of A is due to the fact that the effective field act- 
ing on a particle is larger (than without interactions), 
but not too large to suppress rotations. This finding is 
consistent with the mechanism found in Ref. [Toj . 

Interestingly, we were not able to attribute the in- 
crease in the ratchet effect in systems of dipoles to a 
synchronization phenomenon, i.e., a coupled rotation of 
two neighboring dipolar particles. It is, however, possible 
that such synchronization phenomena occur at thermo- 
dynamic and field parameters that differ from the ones 
investigated here. 

At higher values of A, i.e., stronger dipolar couplings, 
we find a decrease in the ratchet effect. In this region, 
the particles start to aggregate into clusters along the 
direction determined by the constant contribution to the 
external field. As a consequence, the effective field be- 
comes too strong and the dipole moments can follow the 
oscillatory motion of the field less and less, leading to a 
pronounced dampening of any rotations. We note that 
the values of A considered in this work are in the range 
typical for real ferrofiuids (as are the considered densi- 
ties). 

As a somewhat counterintuitive effect, we have found 
that not only the anisotropic dipolar interactions but also 
the isotropic repulsive interactions between the particles 
have a significant influence on the ratchet effect. At con- 
stant dipolar coupling strength, the steepness of these in- 
teractions determines the average distance between the 
particles and thus, the magnitude of the effective local 
field. In this way, short-range interactions can "tune" 
the effective torque. 

In summary our results show that the conservative 
interactions typical of real ferrofiuids strongly influence 
noise- induced phenomena such as the ratchet effect. So 
far, we have not taken into account the fact that the sol- 
vent, which is omnipresent in a ferrofluid, induces addi- 
tional hydrodynamic interactions between the magnetic 
particles. These long-range interactions have been shown 
to play a significant role in translational ratchets (see, 
e.g., Refs. [28|, [29|) and related synchronization phenom- 
ena [30| . The interplay of hydrodynamic and dipolar in- 
teractions in the context of the present ratchet effect will 
be the subject of a future study. 

ACKNOWLEDGMENTS 

We thank A. Engel for the motivation of this work and 
for helpful discussions. Financial support from the DFG 
within the RTG 1558 Nonequilibrium Collective Dynam- 
ics in Condensed Matter and Biological Systems, project 
Bl, is gratefully acknowledged. 

Appendix: Reduced units 

Here, we show how the reduced units used in our study 
are related to the ones used by Engel et al. Q. The 
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latter are denoted by a superscript "f" . For the reduced the time by 
temperature we find 



T* = 



The frequency is related by 



B* 

w 



(A.2) 



(A.l) and the torque by 



(A.3) 



(A.4) 



Using L>t = 0.2, fi* = 1 [cf. Eq. ([TOJ)], and choosing fit = 
B* and = 0.0015, yields T* = 0.2, At* = 0.0001, 
and uj* = 15. This means that T* = T}. 
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